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Abstract 

Exact equations are proposed to describe g-function flows in integrable boundary quan- 
tum field theories which interpolate between different conformal field theories in their 
ultraviolet and infrared limits, extending previous work where purely massive flows 
were treated. The approach is illustrated with flows between the tricritical and critical 
Ising models, but the method is not restricted to these cases and should be of use in 
unravelling general patterns of integrable boundary flows between pairs of conformal 
field theories. 
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1 Introduction 



Since work by A.B. Zamolodchikov more than 20 years ago [1], many examples of two- 
dimensional quantum field theories which flow between different conformal field theories 
in their short and long distance limits have been found. If such theories are placed on 
manifolds with one or more boundaries, then the corresponding boundary conditions must 
also flow, between conformal boundary conditions appropriate to the conformal field the- 
ories sitting at the two limits. To figure out the resulting pattern of combined bulk and 
boundary flows is an interesting problem, with potential relevance to a variety of issues in 
condensed matter physics and string theory. 

Zamolodchikov's original paper (see also [2]) concerned the bulk perturbation of the 
unitary minimal model A4 PtP +i by its 4>\3 operator. For p large, a perturbative calculation 
of the c- function [3] enabled him to show that, for one sign of the coupling, the resulting 
renormalisation group flow interpolates between M pp +i in the ultraviolet and A4 p -i sP in 
the infrared. The generalisation of this approach to the boundary situation is surpris- 
ingly tricky, but has recently been achieved in a paper by Fredenhagen, Gaberdiel and 
Schmidt-Colinet [4], where the g-function or boundary entropy [5,6] was used to identify 
the destination boundary conditions, again for large p. (Even more recently, the same 
pattern of flows has been shown to hold on fluctuating surfaces with boundaries [7].) 

The calculations of [4] are perturbative in 1/p, and they do not give reliable information 
about flows near the bottom of the minimal series. In fact, even at large p the authors of [4] 
had to borrow some non-perturbative information about pure-boundary flows in order to 
obtain a full picture. In the absence of boundaries, Al.B. Zamolodchikov showed how such 
problems could be circumvented in integrable situations through the use of exact, nonper- 
turbative equations of Thermodynamic Bethe Ansatz (TBA) type [8]. These equations 
encode the evolution of a quantity called the effective central charge, c e g, during renor- 
malisation group flows, where c e fr is an off-critical generalisation of the central charge of a 
conformal field theory, which agrees with the c-function used in [1] at fixed points of the 
renormalisation group. The purpose of this paper is to show that a similarly-exact descrip- 
tion of bulk flows with boundaries is possible, at least in cases where the combined bulk 
and boundary theory is integrable. Our starting-point is the exact off-critical g-function 
for massive integrable quantum field theories that was introduced in [9] and further stud- 
ied in [10]. After some background material in section [21 the proposed massless variant 
of the TBA-inspired exact ^-function of [9] is introduced in section [31 together with some 
numerical illustrations of its implications. These results are backed by exact calculations 
of limiting ^-function values in section 0J where we also report some simple perturbative 
checks of our proposal. Finally section [5] contains some conclusions. 

In cases where the bulk remains critical, the use of equations of TBA type to evaluate g- 
functions has a long history, dating back at least to work on the Kondo problem [5] . In this 
respect the main novelty of our result is the demonstration that, for off-critical interpolating 
flows, bulk-induced changes to g-functions can also be accounted for, exactly, through the 
TBA approach. Some motivation for our specific proposal came from a consideration of 
Al.B. Zamolodchikov's staircase model [11] (see also [12-15]). The full set of flows implied 
by this connection is rather rich, and we postpone its discussion to another occasional. We 
have also limited the treatment in this paper to flows between the tricritical Ising and 
Ising models (noting, though, that these cases are of particular interest, being the furthest 
possible from the perturbative limit studied previously). Generalisation to other cases 

1 Though we should mention a previous attempt to use the staircase TBA to study bulk- and boundary- 
induced g- function flows between conformal field theories, reported in [16]. However the equations developed 
there do not fully account for the effects of an off-critical bulk on the ^-functions, and do not agree with 
our results. 
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appears to be straightforward, and we plan to present a more detailed analysis elsewhere. 



2 Background 
2.1 The bulk flow 

The TBA system found in [8] for the tricritical Ising to Ising flow encodes the ground state 
energy E(R) of the interpolating theory on a circle of circumference R via a pair of integral 
equations for two pseudoenergies e\{9) and €2(9). Using the symmetry €2(6*) = e\(—0) these 
boil down to a single equation, for e(9) = e±(9) = €2(—9) : 



e(0) = \ re 9 - [ <j>{6 + 9') L(9') dO' . (2.1) 
2 Jr 



Here L{9) = ln(l + e' 6 ^), <p(6) = 27r C osh(g) > and r = MR with M a parameter with the 
dimensions of mass which sets the (inverse) crossover scale. Then 

E(R) = -^c cS (r) (2.2) 

where 

Cefr(r) = 4r / re e L{6)d6. (2.3) 

The limiting values c e g(0) = 7/10 and c e ff(oo) = 1/2 can be calculated exactly [8], and 
match the central charges of the tricritical Ising and Ising models. Later, the form of 
L(9) in these two limits will be important. As r — > 0, three regions develop where L(6) is 
approximately constant: 

L{9) ~ m(2) = 0.6931 ... for 6 < - ln(l/r) ; (2.4) 

L{0) ~ln((3+V5)/2) = 0.9624... for - ln(l/r) < < ln(l/r) ; (2.5) 
L(9) ~ for 6 > ln(l/r) . (2.6) 

In the opposite limit, r — > 00, there are instead just two regions: 

L{6) ~ ln(2) for 9 < - ln(r) ; (2.7) 

L{6) ~ for » -ln(r) . (2.8) 

These behaviours are illustrated in figure [TJ 
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2.2 The conformal boundary conditions 



With the addition of boundaries, the endpoints of the interpolating flows become boundary 
conformal field theories: the boundary tricritical Ising model in the ultraviolet, and the 
boundary (critical) Ising model in the infrared. The basic (Cardy) boundary states for 
these two models follow from [17], and their physical interpretations were discussed in 
[17] (for Ising) and [18, 19] (for tricritical Ising). For Ising, there are three possibilities, 
corresponding to the boundary spins being fixed up (+), fixed down (— ), or free (/). 
Written in terms of Ishibashi states [20] | 0)), |e)) and | a)) they are [17] 



!(+)> 



V2 



o: 



i 



!(-)) = ^io; 



-~k a)) 



!(/)) = I o»-k» 



(2.9) 



The inner products of these states with the vacuum |0) give the corresponding values of 
the conformal g- function [6] . On Ishibashi states | a )) we have (0 | a )) = 5o a , so 

1 
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I Ising 
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I Ising 



In — = -0.3465 . 



Ising 

We will also treat the superposition of (+) and ( 



V2 
lnl 



(2.10) 



0. 



(2.11) 



-) boundaries, (+)&(—), for which 



hi 



#(+)&(-) Ling ~ ln2 * 



(+) I Ising 



0.3465 . . . 



(2.12) 



For the tricritical Ising model there are instead six options, each labelled, roughly 
speaking, by the value (or values) available to the order parameter (a) at that boundary, 
taken from {—,0, +} [18]. (In the conformal field theory, a becomes the leading spin field, 
with dimensions (3/80,3/80).) These are (-), (0), (+), (-0), (0+), and (-0+), though 
the last of these is traditionally labelled as (d), with d standing for 'degenerate'. The 
corresponding boundary states ared 



!(+)> 


= C 
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»+*/lA»+i7l§» + | 


1) 
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= c 
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1(0)) 


= V2C 
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lo))-^ 1 !^))-^ 1 ! 


!) 


) + v 
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(2.13) 



where 



C 



\J sin(7r/5)/v / 5 



8 8^ 



1/4 



n 



^2 008(^/5) = J(l+VE)/2 . (2.14) 



2 Note, our assignments of the |(+)) and |(— )} states, and of the |(0+)) and |(— 0)) states, are opposite 
to those in [18]. This is to ensure that the one-point function of the spin field a is positive in the presence 
of the (+) boundary, and negative in the presence of the ( — ) boundary, which is more natural, and matches 
the convention adopted in (|2.9|) for the Ising boundary states. We've also corrected a typo in the |(+)) and 
|(— )} boundary states as given in [18]; our states match those given in, for example, [21]. 
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Again adding in the (+)&(—) superposition, the ^-function values we will need, ordered 
by increasing g and expressed in ways that will be useful for comparisons later on, are 



1 , / 1 "1 


-0.6680 . . . 


(2.15) 




-0.3214. . . 


(2.16) 




-0.1868... 


(2.17) 


. . .. + iln2 = 

tricrit 2, 


0.0250 . . . 


(2.18) 


. . .. + iln2 = 

tricrit 2. 


0.1597... 


(2.19) 



^ ^(+) I tricrit ^^(—) I tricrit ^(^0 

ln %)ltricrit = MV2C) = 

ln 5(0+) I tricrit = ln 3(-0)ltricrit = H^C) = 
ln #(+)&(-) Itricrit = m2 3(+)|tricrit = ln %) 
ln %) I tricrit = HV2 V 2 C) = lng {Q+) 

2.3 The boundary flows 

The flows which occur when the two models are perturbed at the boundary alone are 
well-understood. In Ising, the (/) boundary admits a single relevant boundary field (ft 13, 
with dimension 1/2. This breaks the Z2 symmetry of the bulk and can be interpreted as a 
boundary magnetic field. Depending on the sign of the perturbation, a flow is induced to 
the (— ) or to the (+) boundary, as shown in figure [2 

(-) — (/) ► (+) 

Figure [2] Flows from the (/) boundary in the Ising model. 

The (+) and (— ) boundaries have no relevant boundary fields, but one can also consider 
their superposition, (+)&(—). The boundary-condition changing operators correspond to 
013 again, and generate the flow illustrated in figure [3l to the free boundary condition 
[22,23]: 

( + )&(-) ► (/) 

FigureO The flow from the (+)&(—) boundary in the Ising model. 

For the tricritical Ising model the structure is richer [18,19] (see also [24-26]). Including 
the superposition (+)&(—), the full map is shown in figure HI 

^ --(+)&(-) --^ 

, I \ 



1 / 
1 / 



\ 1 
\ 1 



(-)—-(-())-► (0)^-(0+)-m+) 

Figure [4j Boundary flows in the tricritical Ising model. Solid lines (red online) 
show flows induced by a (ftis boundary field. The dashed lines (green online) 
are induced by <fti2 (for the lower two lines) or (ft\\ (for the upper two). Finally, 
the dotted lines (blue online) are induced by a combination of <ft\2 and 1^13. 
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With the bulk conformal, the flows induced by 4>i3, (pi2 and (fin (the solid and dashed 
lines on figured]) are all integrable. However, for integrability to survive when the bulk is 
also perturbed, it is not enough for the bulk and boundary perturbations to be separately 
integrable - they must also be compatible with each other. For the (^13 bulk perturbation 
which leads to the interpolating flow to Ising, this is thought to hold if the boundary 
perturbing operator is also (fi\^ [27], and so it is these combined bulk-and-boundary flows 
that we should aim to treat using the exact ^-function. 



3 Exact g-functions for the interpolating flow 

In [9,10], exact equations were proposed for the off-critical g- function in certain massive 
integrable boundary theories. To give our proposal for the interpolating bulk and boundary 
tricritical Ising to Ising flows, we first set 

, sinh(f + ^) i d , , W/rt -sin(7T2;)/(27r) . , 

(x)(6) = >1 A-r , 4>(x)(0) = =hi(x)(0) = — -77^ — / \ ? 3.1) 

sinh(f-^f) {) 2ird0 v A 1 cosh(6») - cos(vrx) ' v ; 

so that the kernel function (f)(9) in the bulk TBA equation (|2.ip is equal to —(f>(i/2)(9), and 

(fi {x) (9)d9 = -(l-\x\)sgn(x). (3.2) 



Now let e(9) solve the bulk TBA equation (|2.1[) for a system on a cylinder of circumference 
r, and suppose a boundary is placed at the end of that cylinder with a boundary condition 
which depends on a further parameter 9 b - We will propose the following expression for the 
logarithm of a g-function In g(r): 

\ng(r) = In g (r) + In g b (r) (3.3) 

where 

(3.4) 

and 

ln 56 (r) = -Iln(2)+ / (&(0) - </>(20))L(9) d9 (3.5) 
^ ji 

with L(0) = ln(l + e~ e ^) as before, and 

MO) = - <P(i/2)(o - o b ) ■ (3.6) 

(Note, the normalisation of (fib here differs by a factor of 2 from that in [9, 10].) 

The expression ()3.3[) has the same general structure as the exact massive g-function 
introduced in [9], with g\,{r) containing the boundary-condition specific parts of the g- 
function, while g${r) is a 'universal' piece which incorporates the effects of the bulk per- 
turbation on the boundary entropy. However, the new formula involves some significant 
changes too - in particular, the infinite series in (13. 4D contains only odd terms, and all 
rapidity combinations in the kernel functions (fi(9i+9t+i) appear as sums. (This second 
aspect is related to the fact that (|3.4p has been written in terms of the single function 
e{9) = ei(0), rather than e x {9) and e 2 (0) = e x (-9).) 

The infinite series gives an expansion for In go (r) about r = 00, but it converges rapidly 
for all values of r, and can be summed exactly at r = 0, and in various intermediate double- 
scaling limits. Before giving these details, figures [21 [6] and [7] show numerically-obtained 
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plots of In g(r) for b = — 15, and +15. The plots were obtained using the first five terms 
from the series (|3.4p . though truncating to just three terms would have given visually 



indistinguishable results. For all three values of Ob, lng(r) tends to In 5(0+) | tricrit * n ^ e 
far ultraviolet, to Ing^ | Ising i n the f ar infrarecH, and undergoes a transition at lnr ~ 0, 
which is where the bulk crossover occurs. For 6^ = —15, there are two further transitions, 
at lnr ~ ±15, while for 9 h = +15 there is one, at lnr ~ —15. 



0.2 



-0.2- 



-0.4 



-0.6- 



ceff/2 



In(gfree) 



In(g0+) 



ln(g0) 



In(gfixed) 



in(g+) 



-20 



-10 





ln(r) 



10 



20 



Figure[5l The exact g-function flow for 6b — —lb . The flow of c e ff (r)/2, running 
from 0.35 down to 0.25, is also shown, to indicate the location and duration of 
the bulk crossover. Tricritcal Ising g-function values are g0+, g0 and g+; critical 
Ising values are gfree and gfixed. 




-0.6 



Figure [SJ The exact g-function flow for 0^ = . Labelling as for figure O 

respectiveiy. We wiil 



Note, though, that these are equaliy the values of hi; 



'(O-Jltricrit aild ln 5(-)l Mng 

mostly leave this ambiguity implicit in the following, but we will return to it briefly later in this section 
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-20 -10 10 20 

ln(r) 

Figure [71 The exact ^-function flow for 0^ = 15 . Labelling as for figure 



The natural interpretation of the plot for 0& = —15 is that the corresponding renor- 
malisation group flow starts with a pure-boundary transition at lnr ~ —15 from the 
(0+)|tricriticai boundary to the vicinity of the (0)| tr icriticai boundary, with the bulk remain- 
ing near to the tricritical Ising fixed point, then undergoes a bulk-and-boundary transition 
with the bulk flowing from tricritical Ising to Ising while the boundary moves from the 
neighbourhood of (0)| tricritical to the neighbourhood of (/)|i s i ng , before finally making a 
further boundary transition, at lnr rj 15, to (+)| ising- For = 0, there is a single com- 
bined bulk-and-boundary transition, from (0+)|tricrit to (+)|i s ing, at lnr r; 0. The absence 
of an independent boundary transition suggests that this case corresponds to the boundary 
perturbation being zero, and we will give further evidence for this claim in the next section. 
Finally, for Of, = 15 there is a pure-boundary transition at lnr r; —15, from (0+)t r icrit to 
the neighbourhood of (+)tricritj followed by a bulk-and-boundary transition to (+) ising at 
lnr ~ 0, and no further transitions. 



Tricritical Ising: 



Ising: 




Figured] Combined bulk and boundary flows predicted by Q3.3p . (13. 4|) and (13. 5|) . 
Renormalisation group fixed points are labelled by their conformal boundary 
conditions, with the bulk theory for the upper row being the tricritical Ising 
model, while for the lower it is the critical Ising model. The flows marked a, b 
and C correspond to Of, = —15, and 15 respectively. 

These results combine to give the picture sketched in figure El which matches the 
predictions made on the basis of large-p perturbative calculations in [4]. Furthermore, 
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taking the limits 9b — > — oo and 9b — > +00 shows that in addition to the 9b = flow from 
(0+)|tricrit to (+)|isin g , there should be bulk-induced flows from (0) ^critical to (/)|i s i ng , and 
from (+)|tricriticai to (+)|i s ing; as explained in [4], these claims match the results of [28]. 

One caveat, though: as mentioned above, strictly speaking our results cannot dis- 
tinguish between (+) and (— ), nor between (0+) and (—0), as the g-functions do not 
distinguish between these pairs of boundary conditions. Physically it is clear that the 
picture given in figure [8l and its image under a global swap of + for — , must be correct, 
but to resolve the issue within the context of exact ^-function flows alone, one would have 
to track the evolution of the inner products of states other than the ground state with 
the boundary state. We expect that this will be possible using pseudoenergies which solve 
excited-state TBA equations [29,30], but we shall leave the further exploration of this issue 
to future work. 

Finally, we need a proposal for the off-critical deformations of the Z2-symmetric 0i3 
flows which run from (d) up to (+)&(—) and down to (0) in figured! We claim that these 
flows are captured by replacing the formula (|3.5j) for In gb(r) by 

hi 56 (r)= f (0 b (^)-0(2^))ln(l + e- £ W)^. (3.7) 

In other words, we simply add g hi 2 = 0.3465 ... to the logarithm of the previous exact 
g- function. The graphs in figures [6] and [7] are then shifted upwards by this constant, 
and the transitions occur at the same values of r as before, but between a different set 
of conformal boundary conditions, as summarised in figure Again, this matches the 
extrapolation of the predictions of [4] down to p = 4. 




Figure[9] Combined bulk and boundary flows predicted by (13.31) . f|3.4|) and (|3 . T[) . 
The labelling convention for renormalisation group fixed points is as for figure[8l 
The flows marked d, e and f correspond to 9b = — 15, and 15 respectively. 



4 Exact and numerical tests of the proposal 
4.1 Exact limiting values of the g-function 

We first deal with the universal factor lng (r) defined by equation (|3.4p . From (|3.2p . 
f R (j>(9)d9 = 1/2, and so 

I d9 1 ...de 2j ^ 1 ^(9 1 + 9 2 )H02 + 9 3 )...^(9 2 ^ 1 + 9 1 ) = ^- (4.1) 

with the product of the kernel functions <p tending to zero exponentially outside a region 
of order 1 about the origin. In the infrared, only this latter property is needed: from (|2.8p . 
each factor l/(l+e(0j)) — ► for Oi 3> — ln(r). In particular this holds in the neighbourhood 
of the origin where the product of the kernel functions is significantly different from zero. 
Hence all terms in the series (I3.4p tend to zero as r — > 00, and g {r) — > 0. 
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In the ultraviolet, via (|2.5p . the pseudo-energies tend to constants in the central region 
-ln(l/r) < « ln(l/r), with e" e W -»■ (l+y/E)/2, l/(l+e e W) -» (V5-l)/2. Combined 
with (|4.ip . this implies 



_ x 1 + e ^) • • • i + e ^-0 + + *») ■ ■ ■ + ^ 2 X2J " (42) 

where x = (v5— 1)/4, and so 

1^ 1 9 , , 1 /l+x\ 1 / 2 \ 
lim InflbW = 2 E 2" I ^ = 4 ln ( T^J = 4 ln + 7EJ = ' 15972912974 

(4-3) 

The boundary condition dependent piece can be treated by rewriting (|3.5p so as to split 
In g b (r) into three parts: a constant lng bl , a parameter- independent piece In g b2 (r), and a 
^-dependent piece ln g b3 (r, Of,): 

ln 5 6 (r) = ln g bl + ln 562 (r ) + ln c/ 63 (r, 6> 6 ) (4.4) 

where 

ln 561 = -^ln2, (4.5) 



\ng b2 (r) = / fa (3/4 )(0) - <I>(20))L(0) d8 , (4.6) 
ln<7 fe3 (r,0 fe ) = -[ <p (1/2) (8-8 b )L(9)d8. (4.7) 



Note also, from (|3.2p . that 



/ 



(0(3/4)^-0(2^)^ = -^ (4.8) 
0(i/2)(^)^ = ^. (4.9) 



The decay properties of 4>( x ) (&) mean that the support for the integral (|4.8p is concentrated 
near to 8 = 0, while that for (|4.9p is concentrated near to 9 = 8 b . Combined with the 
asymptotic behaviours of L(0) recorded in equations (|2.4p - (|2.8p . these results allow the 
various limiting values of lng(r) to be computed. 

1. In the far infrared limit {r — > 00, 8b fixed}, lng (r) — ► 0, and L{8) —* in the regions 
where the integrals (14.61) and (14. 7p can receive contributions, so 

ln<7(r) - ln 561 =-|hi2 = In <7 {+) | Ising • (4.10) 

2. In the far ultraviolet limit {r — ► 0, 8b fixed}, L{6) acquires a constant value in the 
whole region where the integrands in (|4.6p and (|4.7p are significantly different from zero. 
The integrals (|4.6|) and (j4.7H therefore cancel in the limit, and 

\ng(r) -+ ln ffo (0) + ln 5w = | ln(l + ^) - ± ln2 = | ln(± + = ln %+) ^ . (4.11) 

3. If 6 < and 6 < lnr < 0, then for 8 « 6 , L(0) « hi 2 from ([23]), while for 9 « 0, 
L(0) ps In((3+V5)/2) from (|23|) . Hence ln^ ±ln2 and ln 5fe2 « -i ln((3+V5)/2), and 

ln 5 (r) - |ln(l+^) - \ In 2 - \ In((3+V5)/2) + |ln2 = ln 5(0) | tricrit . (4.12) 
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4. If 9 b < and < lnr < -9 b , then for w 9 b , L{9) ps In 2 from ([12]), while for fl»0, 

~ from (|2.8p . Hence mg fe3 ~ | In 2 and hi</ fe2 ~ 0, and 

ln 5 (r) -» -i In 2 + ± In 2 = = ln 5(/) | Ising ■ (4.13) 

5. If 6 > and -9 b < lnr < 0, then for ^0, L(6») ps ln((3+\/5)/2) from (|23]) . while 
for 6» ps 6 , L(0) ps from <|2^j> . Hence ln# b2 ps -| ln((3+\/5)/2) and ln# 63 ps 0, and 

ln 5 (r) -> |ln(l+^) - ± In 2 - \ ln((3+^)/2) = ln 5(+) | tricrit . (4.14) 

To make the statements of 3, 4 and 5 precise, they should be considered as double- 
scaling limits: for example, for 3 one could fix two constants 9 b and f with 9 b < lnr < 0, 
and set r = f p , 9 b = p9 b ; then (|4.12p holds in the limit p — ► oo. If instead is kept fixed 
and r is varied from to oo, then the case 9 b <C 0, figure[5l is covered by 1, 3, 4, 2; the case 
9 ss 0, figure O by 1, 2; and the case 9 b 3> 0, figure by 1, 5, 2. We have thus confirmed 
analytically the previously-observed numerical results, and justified that our conjectured 
equations are indeed consistent with the flow patterns depicted in figures [8] and [9j 



4.2 Comparisons with conformal perturbation theory 

The bulk perturbation which induces the flow from the tricritical Ising model to the Ising 
model corresponds to the addition of a term A f 4>is(x, x) d 2 x to the action of the tricritical 
model, where the bulk coupling A has dimension (mass) 4 / 5 . If the unperturbed conformal 
boundary condition (a) supports the boundary field 4>\z{x), the addition of a boundary 
perturbation p J 4>i3(x) dx can also be considered, where [i is the boundary coupling, with 
dimension (mass) 2 / 5 . (For the tricritical Ising conformal boundary conditions featured on 
figure H (-0), (d) and (0+) do support this field, while (-), (0), (+) and (+)&(-l) do 
not.) A g- function as evaluated in conformal perturbation theory should therefore have 
the expansion 

oo 

\ug{\^,R)= c^ n (pR 2 ^r^R 4/5 ) n ■ (4-15) 

m,n=0 

In general this is a regular series in powers of i? 2//5 , reducing to a series in when 
p, = 0. At large R, the function defined by (|4.15p will typically develop a linear behaviour, 
with ]nQ(h, A, R) ~ —fMR where / is a free energy per unit length, which we choose 
to measure in units of the inverse crossover scale M. Our exact equations, by contrast, 
yield 'subtracted' (^-functions from which this term is absent in the infrared, and instead 
reappears as an irregular term in the ultraviolet [32,33]. They are also expressed in terms 
of M and the boundary parameter 9 b , rather than A and p. The relation between A and 
M is known [8,31]: 



A = kM 4 / 5 , k = ^ r\l S^ Q \ = 0.14869551611... (4.16) 

2 v / 2(37r)V5Y r(3/10) 1 7 

and on dimensional grounds it must be possible to write p as 

)M 2 / 5 (4.17) 

where v is some dimensionless function of 9 b . The (/-function defined by (|3.3|) should thus 
have the following expansion about r = MR = : 

oo 

In g(r, 9 b ) = In g (r) + In g b (r) = fr + £ c$ n (ur 2 / 5 ) m ( K r 4 / 5 ) n (4.18) 

m,n=0 
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where Cq"J is equal to the logarithm of the conformal ^-function gi a \ for the (a) boundary 
condition. For all values of 9 b we expect 

eg = (4.19) 

since c[q is proportional to the one-point function of the perturbing boundary operator 
on a disk with the vacuum field at its centre [33] , and vanishes in a unitary theory such as 
this one [9@. Previous examples suggest that In g (r) will not contribute to the irregular 
term fr [9, 10,33] ; assuming that this holds true here too, the value of / can be calculated 
from (|3.5p as in [32,33], with the result 

m) = h' eb - its • ( 4 - 2 °) 

Finally, the first bulk-induced coefficient in the expansion of In g(r) is [10] 

(a) _ -6(1-3^,3^/2) ((ft | (Op) 



2(2tt) 1 -^ (0|(a)) 



(4.21) 



where B{x,y) = Y(x)Y(y)/Y(x+y) is Euler's beta function, eft is the bulk perturbing field, 
<fti3 in this case, and is its scaling dimension, here equal to 6/5. The inner products 
(<ft |(a)) and (0|(a)) = g,-, can be read from (|2.13p . bearing in mind that the Ishibashi 
states in those formulae have been labelled by the conformal dimensions of their Virasoro 
representations, which are half the scaling dimensions of the corresponding bulk fields. 

The bulk TBA equation (|2,ip was solved numerically for 101 evenly-spaced values of 
r 4/5 runnm g f r om 0.0005 to 0.1255, discretising the 9 axis to 1520 points between 9 = —50 
and 9 = 50 and using extended (20 decimal digit) precision in GNU Fortran 95. The 
resulting estimates for the pseudoenergy e(9) were then used to compute lng(r) via (|3.3p . 
summing the series (I3.4p for In go(r) to 12 terms, and evaluating the ^-dependent part 
lngb(r) from (|3.5f) for values of 9b ranging from —2.5 and 2.5. (Were accurate results to 
be required for a larger range of 9b, care would have to be taken to decrease the values of 
r 4//5 used for the fits, to avoid their contamination by the intermediate plateau values of 
\ng{r) which appear as \9b\ increases, as on figures [5] and [71) 

As a first check of our numerics, we made a least-squares fit of the function c e s(r) 
defined by f|2.3|) to a regular expansion in powers of r 4 / 5 plus a single 'antibulk' term 
proportional to r 2 , finding coefficients which matched those reported in [8] to the full 
accuracy claimed there. 

Then, for each value of 9b, the numerically-obtained In g{r, 9b) was fitted to a series in 
r 2 / 5 plus a single term proportional to r, as in (|4. 18[) : 

oo 

In g(r, 9 h ) = ]T d k {9 b ) r 2k ' b + e(9 b ) r . (4.22) 

fc=0 

If the match with conformal perturbation theory is to hold, we should have 

[k/2] 



d k (9 b ) = Yl ci k^ k ~ 21 ^ > <°b) = f(9 b ) • (4.23) 

1=0 



The constant term c?o(#b) obtained from the fits matched the value predicted by (|4.1ip . 
namely ln<^ 0+ ^ | tricrit , to at least 10 digits for the whole range of 9b- Furthermore, d\{9b) 



4 Note, if c 10 did not vanish, then the ^-theorem, which states that g decreases for all pure-boundary 
flows in unitary models [6,34], would be violated for one or other sign of the boundary coupling. Conversely, 
the non-vanishing of co,i [10] is an easy way to see that the g-theorem can be violated when the bulk flows, 
even in a unitary theory (see also [35]). 
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was zero to the same accuracy, in line with (|4.19p . In figure PTOa. the values of e(9b) found 
from our fits are compared with the exact predictions from (|4.20|) ; the good agreement 
supports our claim that In (70(f) does not contribute to this linear term. 



5H 

2\ 
F 



4- 

d 3 2 ~ 



-2 



-2 

-4- 



Figure [TOk : A comparison of the exact pre- 
diction ([420]) for f(9 b ) (dotted line) with val- 
ues of e(9 b ) from the fit (|4~22ll (points). 



Figure [TUb: Estimates of d^Ob) from the fit 
(gjgg) . The estimate at 9 b = is 5.24xl(T 8 , 
consistent with to our numerical accuracy. 



Next, in figure [TUb. we show the values of d^{9b)- The apparent zero of this function at 
9b = suggests that this point should correspond to \i = in (|4.18|> . where all odd terms 
in the regular series should vanish. This is consistent with the fit at 9b = 0, which is 



In g(r, 0) 



-0.1868444605395000 + 0.1464465945456 r - 8.429x 10~ 12 r 2/5 
- 0.2038867755577 r 4 / 5 + 5.24x 10~ 8 r 6 / 5 - 0.008541571 r 8 / 5 



+ 3.68xl0~ 6 r 2 



0.00209 r 



12/5 



+ ... 



(4.24) 



Supposing that the odd terms are exactly zero for 9b = 0, a more-constrained fit to a 
regular series in powers of r 4 / 5 plus a term linear in r gives the result 

lng(r,0) = -0.1868444605395363 + 0.1464466094005 r - 0.2038867770734 r 4/5 

- 0.008541178 r 8 / 5 - 0.0020624 r 12 / 5 + 0.00151 r 16 / 5 - 0.0004 r 4 + . . . (4.25) 

For this case, c?2 is known exactly, since with /j and hence v = 0, cfo = Cg^ft and can be 
calculated from (|4.21j) and (|4.16j) . For the (0+) boundary with 9b = we thus have the 
exact predictions 



do 



In 



9(o+) 



-0.1868444605395326. 



di = 



£(-1/5,3/5) /« 
2(2vr)-V5 r? 3 



-0.2038867770751855 . . . 



«3 

e 







/(0) = 0.1464466094067... 



(4.26) 
(4.27) 
(4.28) 

(4.29) 
(4.30) 



all of which are reproduced very well by the fits (|4.24|) and (|4.25p . 

The limits 9b — > — 00 and 9b — > +00 admit similarly-simple checks. A consideration 
of figure [8] and the results from subsection 14.11 shows that if these limits are taken first, 
keeping r finite, and r is only then allowed to vary, then the resulting equations should 
describe the bulk-induced flows (0)| tr icrit - * (/)|ising and (+)|tricrit - * (+)ising respectively. 
Neither of the UV boundary conditions for these flows admit a ^13 boundary field, so the 
logarithms of their g-functions should have regular expansions about r = in powers of 
r 4//5 , with a coefficient di of r 4//5 that can be predicted from (I4.2ip . 
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In the limit 6^ — > — oo, making use of (|4.9p and (|2.7|) . equation (|3.5p reduces to 

ln^(r)|, i;= _ 00 = / {<l, im {6)-<f>{7B))L{0)dB (4.31) 



while equation (|3.4p for lng (r) is unchanged. With the numerical work as before, the fit 
to a series in powers of r 2 / 5 , together with a linear term, was 

^9(r)\e b =-oc = -0.3214826953191443 - 0.3535533993721 r - 5.178xl0~ 12 r 2/5 
+ 0.5337825131495 r 4/5 + 3.04x 10~ 8 r 6/5 - 0.01741761 r 8/5 

+ 1.98x 1(T 6 r 2 + 0.0132 r 12/5 + . . . (4.32) 

and the more-constrained fit to a regular series in powers of r 4 / 5 plus a linear term gave 

ln^r)!^.^ = -0.3214826953191671 - 0.3535533905994 r + 0.5337825122412 r 4/5 

- 0.017417394 r 8/5 + 0.0133024 r 12/5 - 0.00130 r 16/5 - 0.0008 r 4 + . . . (4.33) 

These results can be compared with the exact predictions for the first few coefficients for 
the bulk-induced flow from the (0) boundary: 

d = In g (Q) = -0.3214826953191634 . . . (4.34) 

&x = (4.35) 
B(-1/5,3/5)kt7 

d 2 = v / ' ' ' — '- = 0.5337825122395085 . . . 4.36 

2(2tt)- 1 /5 v ; 

d 3 = (4.37) 
e = _ 1 = -0.3535533905932 . . . (4.38) 

Again, the agreement is very good. It is also straightforward to check analytically that this 
g-function interpolates between the desired values - the argument is essentially covered by 
cases 3 and 4 of the last subsection. 

For the (+)|tricrit — ► (+)ising flow expected to arise in the Q\, — > +oo limit the story is 
very similar. Since, by (|2.8p . L(9) — > as 9 — > +oo, there is this time no modification to 
the constant term in (13. 5p . which becomes 

ln 56 (r)| 05=+oo = -iln(2)+ / (<f> m) (9) - ^(2fl))L(6») dQ (4.39) 

1 JR 

with the expression for In g Q (r) again unchanged. The fits for the expansion coefficients 
(apart from the constant term) are therefore the same as before, and it is straightforward 
to check that these match expectations from perturbation theory for this situation. The 
same also holds for the set of ^-symmetric flows predicted by our second proposal, (|3.7p . 
and so we will leave the details to the reader. 



5 Conclusions 

In this paper we have proposed an extension of the exact off-critical g-function equations 
of [9,10] to cover situations where the bulk field theory retains massless degrees of freedom 
even in the far infrared, and therefore interpolates between two different conformal field 
theories. While our proposals are still conjectural, they have passed a number of non-trivial 
checks against perturbation theory, leaving us in little doubt that they are correct. Never- 
theless, a first-principles derivation from field-theoretic considerations would be valuable, 
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as would an understanding via lattice models, as obtained for the TBA equations for the 
bulk interpolating flow of c e ff(r) in [36]. However this may be a hard task, and indeed 
previous attempts to derive exact equations for exact (^-function flows have not been en- 
tirely successful [370. It would also be interesting to see whether a process of analytic 
continuation could be used to relate massive and massless (^-functions, though again this 
might be delicate, given that the functions are initially defined via a sub-leading term in 
the asymptotic behaviour of cylinder partition functions. 

One feature of our results is the exact equality, up to a constant factor, of the g- 
functions for various a-priori different flows. It seems likely that this can be understood 
through an extension of the defect-related work of [38] to theories off-critical in the bullfl 
and it would be interesting to explore this further. 

The tricritical to critical Ising flow is interesting in its own right (see for example [8,39]), 
but the main reason for concentrating on this particular case in this paper has been its 
relative simplicity, which has allowed us to make a detailed check of the feasibility of our 
approach, and to illustrate the main ideas without too many distracting complications. As 
mentioned in the introduction, we expect that the general method will be of much wider 
applicability, and we hope to return to its further applications in the future. 
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